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Abstract. We study a gauge invariant order parameter for deconfinement and the chiral condensate in 
SU(2) and SU(3) Yang-Mills theory in the vicinity of the deconfinement phase transition using the Landau 
gauge quark and gluon propagators. We determine the gluon propagator from lattice calculations and the 
quark propagator from its Dyson- Schwinger equation, using the gluon propagator as input. The critical 
temperature and a deconfinement order parameter are extracted from the gluon propagator and from 
the dependency of the quark propagator on the temporal boundary conditions. The chiral transition is 
determined using the quark condensate as order parameter. We investigate whether and how a difference 
in the chiral and deconfinement transition between SU(2) and SU(3) is manifest. 
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1 Introduction 

Two of the most characteristic features of QCD are at the 
same time two of the most elusive ones: Confinement and 
chiral symmetry breaking. Of particular interest is the de- 
pendency of these phenomena on temperature and density. 
For vanishing current quark masses chiral symmetry gets 
restored by a phase transition above a critical tempera- 
ture. On the other hand in the limit of infinitely heavy 
quarks a phase transition from a confining phase to a de- 
confined phase takes place. Of course, in full QCD with 
physical quarks none of the aforementioned limits is ap- 
propriate. Indeed, at zero chemical potential the would-be 
order parameters for both transitions show a rapid change 
signaling a well-marked cross-over region [T] . Furthermore 
it seems remarkable that the temperature ranges for both 
cross-overs are notable close to each other [2j[3l|4] . This is 
in contrast to e.g. the case of adjoint quarks where both 
temperatures differ by a factor of almost eight [5]- 

Thus, it warrants to investigate the chiral and de- 
confinement transition for fundamental quarks closely, to 
eventually understand both properties. To avoid the prob- 
lems involved with dynamical quarks, a useful laboratory 
is the quenched approximation where only valence quarks 
are present. In this case the deconfinement transition can 
be cleanly defined by usage of the Polyakov-loop as order 
parameter, and both the chiral restoration and deconfine- 
ment temperature coincide within available accuracy [6]. 



Useful tools to address these questions are the so-called 
dual observables [7t[51[^fTUirTTirT^rT^[T^P^l[TBiri7] . These 
are obtained by modifying the temporal boundary con- 
ditions of the valence (test) quarks, without altering the 
dynamics of the theorjQ. Dual observables are sensitive 
to the spectral properties of the Dirac operator, and thus 
do encode both the confinement and the chiral properties 
of quarks. In particular, they are not only able to dis- 
tinguish the low-temperature from the high-temperature 
phase, but also phases where chiral symmetry is broken 
but confinement is no longer present |12j . 

Dual observables have been first defined using lattice 
gauge theory [8] , but have turned out to be quite expensive 
to determine. In particular, the large distances character- 
istic of confining physics and the small masses relevant to 
chiral symmetry are hard to reach. On the other hand, it 
is also possible to determine dual quantities using func- 
tional methods |131ll7j . In this case, there is no limit nei- 
ther on distance nor on the quark masses. The challenge 
is, however, that truncations are necessary. In this work 
we combine the best of both worlds by a combination of 
lattice and continuum methods [TSIfbUfrg] . 

The goal of our study is twofold. On the one hand 
we investigate (potentially gauge-dependent) mechanisms 
linking deconfinement and chiral symmetry restoration. 
The study presented in this work in particular highlights 
the role of the longitudinal, electric part of the Landau 
gauge gluon propagator in this respect. On the other hand 
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Such an alteration would correspond to the introduction of 
an imaginary chemical potential [18] as discussed in Ref. [17] . 
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our study serves as an important intermediate step to- 
wards an analysis of the QCD phase diagram at non- 
vanishing chemical potential. Due to the notorious sign 
problem of lattice QCD in that realm there is great de- 
mand for other methods beyond model calculations. Func- 
tional methods like Dyson-Schwinger equations [20[[2"T] and 
the functional renormalization group [22] are ready to fill 
this gap. To provide for meaningful results at non-zero 
chemical potential it is, however, first necessary to gain as 
much insight as possible into the fidelity of truncations at 
zero chemical potential, where one can compare with the 
results from lattice QCD. As a result of this work we find a 
clear (gauge-invariant) signal for the deconfinement tran- 
sition from the dressed Polyakov-loop determined from the 
Dyson-Schwinger equations in agreement with the infor- 
mation we extract from the lattice gluon propagator. 

This paper is organised as follows: Our method to 
determine the magnetic and electric part of the Landau 
gauge gluon propagator on the lattice are discussed in sec- 
tion [2] Though the gluon propagator is gauge-dependent, 
it is found that it can be used to obtain gauge-invariant 
properties, like the critical temperature. The gluon prop- 
agator is also one of the key ingredients into the Dyson- 
Schwinger equation for the quark propagator. In princi- 
ple the gluon propagator could be obtained using Dyson- 
Schwinger equations il9,23,24ll2"5]. but this turns out to 
be a formidable task [19]. Instead we use our result from 
lattice gauge theory as input, as described in section [3] 
From the Landau gauge quark propagator we then extract 
the gauge- invariant order parameters for the chiral and de- 
confinement transition. This procedure has been used pre- 
viously in [T3"lll4| . However, the then available lattice data 
(T9"] for the gluon propagator have been very coarse on the 
temperature axis, and only available for SU(2). We expand 
here to a much finer grid in the temperature domain, while 
at the same time employing much larger physical volumes 
on the lattice. Preliminary results of these calculations can 
be found in [26]. In addition, we study the physical case 
of SU(30- The results are summarized in section 2] 



2 The gluon propagator from lattice 
calculations 



2.1 Simulation details 

To determine the gluon propagator between zero and roughly 
two times the critical temperature (T c — 277 McV for 
SU(2) and T c = 303 MeV for SU(3)), lattice gauge theory 
implementing the Wilson action is used, details can be 
found in Qj2[28| . Details of the simulation parameters can 
be found in appendix [A] in particular its table [T] There 
also systematic errors will be discussed. This investigation 
extends previous works in the same temperature range [191 
|2"9"] . and at higher temperatures [50] . 



At finite temperature, the Landau gauge gluon propa- 
gator D^ u is described by two independent dressing func- 
tions, Dt and Dl, 

DfM = Pl{ P )D^ b {Lol,p 2 ) + P^Df (o;> 2 ) 

PJM = (i - <W(i - M (<W - Vj §f) 

where P^ v and P^ v are projectors transverse and longitu- 
dinal to the heat bath. The color-dependency of the prop- 
agator has been evaluated in lattice gauge theory, and has 
always been found to be consistent with a color-diagonal 
propagator [28] . Thus there remain two scalar functions 
Dt = Zt/p 2 (transverse to the heat bath's four- velocity, 
chromomagnctic) and Dl = Z^/p 2 (parallel to the heat 
bath's four- velocity, chromoelectric) , with Zt = Zl at 
zero temperature, for the description of the gluon. Here p 
denotes always the four- momentum with p 2 = ui 2 + p , 
which coincides with the three-momentum for the soft 
modes uj p = 0. Only this soft mode will be discussed in 
detail here. The hard modes are discussed separately in 
appendix [B[ and found to be very well approximated by 
D(u; 2 p> p 2 ) = D(0,u 2 +p 2 ). 

The gluon propagator, as the gluon itself, is only well- 
defined in a fixed gauge. For the purpose of calculating 
gauge-invariant quantities the choice of gauge is irrele- 
vant, but since gauge-dependent intermediate results will 
be transferred from one method to another here, a non- 
perturbatively unambiguous definition of the gauge seems 
necessary. We choose here the minimal or average- B Lan- 
dau gauge [3TU32] . For the discussions of other choices see 
[3"2Tl3"3"ll3"4Tl3l)] . However, the gluon propagator at the vol- 
umes employed here is depending only marginally on this 
choice [32 , 34J . In particular, the effect is likely only of the 
same order as other systematic effects of the present cal- 
culations. Therefore, it can also be used for similar other 
calculations, independent of the particular choice between 
all non-perturbative realizations of the Landau gauge. The 
methods used to fix this gauge at finite temperature can 
be found in [TH] for SU(2), and the extension used to ad- 
dress SU(3) in [3H]. In [TW2"8"] also the methods used to 
determine the gluon propagator are presented. Note that 
the method to determine the value of the gluon fields em- 
ployed in general and in particular in [19[ is only applicable 
to configurations with positive real part of the Polyakov 
loop [2S]. Therefore, only such configurations have to be 
included in the Monte-Carlo average which satisfy this 
condition. This is performed by filtering at all tempera- 
tures for this condition. 

In addition to the gluon we determine the ghost prop- 
agator Dq, related to its dressing function by 



D G (p) 



G( P ) 



(1) 



2 First explorative studies of the gluon propagator in the 
unquenched case can be found in [27] . 



and found previously to be color-diagonal [28] . In contrast 
to the gluon propagator, the ghost is expected to depend 
significantly on the gauge choice for momenta p < 1 GcV 
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SU(2) transverse gluon propagator below T 



SU(3) transverse gluon propagator below T 
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Fig. 1. The (soft mode of the) transverse part Dt of the gluon propagator at temperatures below (top panel) and above (middle 
panel) the phase transition. The bottom panel shows both results together. On the left results for SU(2) and on the right results 
for SU(3) are shown. Lines are drawn to guide the eye. Momenta here and hereafter are aligned along the a;-axis. 



32,33,34]. Hero, the only reason for its investigation is to 
study the prediction of functional methods [T^] that the 
gauge-fixing dynamics, which is essentially given by the 
ghost propagator, is insensitive to the phase transition. 
This is indeed the case as already indicated by the results 



of Ref. [TO] and further confirmed by our results for the 
ghost discussed in appendix [U] 

All results presented arc renormalizcd such that at 2 
GeV the dressing functions in all cases equal 1, for con- 
venience. The renormalization factor is determined by the 
average value of the propagators at the closest momen- 
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SU(2) longitudinal gluon propagator below T 



SU(3) longitudinal gluon propagator below T 
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Fig. 2. The (soft mode of the) longitudinal part Dl of the gluon propagator at temperatures below (top panel) and above 
(middle panel) the phase transition. The bottom panel shows both results together. On the left results for SU(2) and on the 
right results for SU(3) are shown. Lines are drawn to guide the eye. 



turn values above and below 2 GeV in the x-y-plane. This 
procedure entails that the renormalization constants are 
temperature-dependent. They vary by about 10% for the 
transverse gluon and the ghost, but by up to 50% for the 
longitudinal gluon in the investigated temperature and 
discretization range. 



2.2 Results 

The finite-temperature results for the transverse part Dt 
of the gluon propagator are shown in figure [TJ In the 
left panel we display the results for SU(2) and in the 
right panel corresponding results for SU(3). The temper- 
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| Electric screening mass for SU(2) | | Electric screening mass for SU(3) | 
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Fig. 3. The electric screening mass as a function of temperature (top panel), and zoomed in the phase transition region (middle) 
panel. The corresponding susceptibility %t — <9-Dz,(0)~ 1//2 /<9T is shown in the bottom panel. On the left results for SU(2) and 
on the right results for SU(3) are shown. Note the respective different range for SU(2) and SU(3) for the susceptibility. Lines 
are drawn to guide the eye. 



ature behaviour of each momentum mode measured can 
be read off the 3d-plot at the bottom of each panel. In 
general we observe a relatively smooth variation of all 
large momentum modes of the propagator across the phase 
transition in agreement with previous expectations [19]. 



These modes behave similar for SU(2) and SU(3). Only 
in the infrared some more drastic variations can be seen. 
These appear to be systematic close to the phase transi- 
tion temperature and are somewhat more pronounced for 
SU(3). Although these variations are clearly not statisti- 
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cal in nature, it is not clear whether they are genuine. 
In particular, they may be volume effects, since increased 
long-wavelength thermal fluctuations close to the phase 
transition arc not faithfully represented by the volumes 
presently employed. This possibility is supported by the 
fact that the effect is stronger in the case of SU(3) than 
of SU(2) where the employed volumes are smaller. In any 
case, more refined and systematic studies are necessary to 
explore the nature of these variations. 

In general we observe that the chromomagnetic, trans- 
verse sector of Yang-Mills theory shows no pronounced 
reaction to the phase transition. The only essential effect 
is that with increasing temperature the propagator be- 
comes stronger infrared suppressed and develops even a 
maximum at p sa 0.5 GeV for the largest temperatures 
investigated here. This is indeed expected given that the 
propagator in the infinite-temperature limit of the dimen- 
sionally reduced theory in this gauge exhibits a clear max- 
imum |39| . 

The situation is rather different for the chromoclectric, 
longitudinal part Dl of the gluon propagator, shown in 
figure [5] From the plots it is immediately visible that the 
longitudinal gluon reacts strongly to the phase transition, 
in accordance with previous observations [19] . For temper- 
atures below the phase transition we observe a dramatic 
increase of the infrared part of the propagator close to the 
phase transition which leads to an even more dramatic 
decrease shortly above T c . Also this dramatic variation 
seems to be more pronounced in the SU(3) data. This be- 
haviour is most easily visible from the electric screening 
mass 

m L =D L {p^0r 1 '\ 
and its associated susceptibility 

_ dm L 
XT - dT , 

which are shown in figure [3] This mass drops when moving 
towards the phase transition, reaches a minimum at the 
phase transition, and then quickly increases. In fact, the 
sensitivity is sufficiently strong that the gauge-invariant 
critical temperature can be determined independently from 
the chromoelectric screening mass just by extracting it 
from the plot. Indeed, this temperature coincides with the 
independently determined one using the Polyakov loop as 
an order parameter [37U38j within the resolution of the 
temperature grid employed here. 

In the middle panel of figure [3] we zoom into the tem- 
perature region around the phase transition. From the 
available data, it appears that the screening mass in the 
SU(3) case is changing more rapidly across T c than for the 
SU(2) case. It may be that this more rapid change signi- 
fies the first-order nature of the SU(3) transition, whereas 
in the SU(2) case the behavior is smoother as may be 
expected for a second-order phase transition, despite the 
larger physical volumes. This interpretation may also be 
supported by the much larger susceptibility in the SU(3) 
case, shown in the lower panel of figure [3] However, at 
present this interpretation may only indicate a possible 



scenario and certainly has to be checked in a more de- 
tailed analysis [40] . In addition, a careful volume study is 
needed to check for the presence of long range correlations 
in the form of a vanishing electric screening mass of the 
SU(2) gluon at T c . 

Together, our results for the chromoelectric and chro- 
momagnetic part of the gluon propagator indicate that the 
dominant response to the phase transition occurs in the 
chromoelectric sector. This finding can be related to the 
behaviour of certain dimension- two condensates [4T1I42) . 
These condensates arc effectively given by the difference 
between the traces of the transverse and the longitudi- 
nal gluon propagator [33]. Thus the asymmetry observed 
in these condensates at finite temperatures [JT] can be 
attributed essentially to the behavior of the longitudinal 
gluon propagator alone. 

These comments complete our investigation of the tem- 
perature dependent gluon propagator. In the next section 
we will see that the gluon data are an important input 
into the Dyson-Schwinger equation for the quark propa- 
gator and therefore serve as a basis for our investigation 
of the order parameters for the chiral and deconfinement 
transition of quenched QCD. 

3 Order parameters extracted from the quark 
propagator 

3.1 The conventional and the dual quark condensate 

The temperature dependent general expression for the in- 
verse quark propagator S -1 is given by 

S^ 1 (p 1 u> p ) = ij 4 uj p C(p,ujp) 

+i-f i p l A(p,uj p ) + B(p,uj p ) , (2) 

with vector and scalar quark dressing functions C, A, B. A 
further tensor component proportional to <r^ v is possible 
in principle but is negligible in the presence of a pure vec- 
torial quark-gluon vertex |21j . as used here. The momen- 
tum arguments arc given in terms of the three momenta 
p and generalised Matsubara frequencies uj p 

LO p (n t , ip) = (2nT)(n t + (p/2-ir). (3) 

These correspond to t/(l) valued boundary conditions for 
the quark fields, i.e. ip(l/T,x) = e lv ip(0,x). The angle 
ip varies between (p G [0,27r[, with (p = for periodic 
and <p — 7r for the usual, physical antiperiodic boundary 
conditions for the fermionic quarks. 

Given the momentum behaviour of the non-perturbative 
dressing functions C, A, B one can extract a (^-dependent 
quark condensate from the propagator according to 

= Z 2 N c tY, J tSj tr c S(p, cjM) . (4) 

The conventional quark condensate is obtained for the spe- 
cial case ip — 7r and multiplication of this expression with 
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Z m . In the limit of vanishing bare quark masses it is an 
order parameter for the chiral phase transition. 

The corresponding dual observable, the dual quark 
condensate or dressed Polyakov loop, is obtained by a 
Fourier-transform of the (^-dependent condensate with re- 
spect to the winding number n, 



2tt 



(5) 



and specialising to the case n = 1. The dressed Polyakov 
loop Si is sensitive to the breaking and restoration of 
center symmetry [DIE] and therefore serves as an order 
parameter for the deconfinement transition. 

These two order parameters, Eq. <^ with tp = n and 
Eq. (|5|), can be extracted from the dressed, temperature 
dependent quark propagator calculated from functional 
methods as detailed in Refs. [T51[T4l[T7] . In the following 
we will update the calculation presented in Ref. [14] and 
in addition also extend the studies to the case of the gauge 
group SU(3). We will investigate whether and how the bet- 
ter temperature resolution of the lattice data for the gluon 
propagator discussed in the previous section leads to im- 
proved results for the order parameters. We also compare 
the results obtained for gauge group SU(2) and SU(3). 



3.2 Truncation scheme for the Dyson-Schwinger 
equations for the quark propagator 

The Dyson-Schwinger equation for the quark propagator 
at finite temperature T is given by 



S- 1 (p, Up ) = Z 2 S^ip, uj p )~C f ^ g 2 T 
X7 M S{k,L> k )r v {k,L>k,p, Up) 



(2^)3 



(6) 



Here Z? M „ denotes the (transverse) gluon propagator in 
Landau gauge and we have introduced a reduced quark- 
gluon vertex _T„, by defining = ig^F v . The bare 

quark propagator is given by Sq (joi) — i-f ■ p+ Z m m(fJ, 2 ), 
where m(fi 2 ) is the renormalized current quark mass. The 
wave function and quark mass renormalization factors, Z 2 
and Z m , are determined in the renormalization process. 
The ghost renormalization factor Z% is canceled by a cor- 
responding factor in our model for the quark-gluon ver- 
tex discussed below. Furthermore we used Z\ = 1 for the 
renormalization factor of the Landau gauge ghost-gluon 
vertex. The quark dressing functions A(p, u) p ), B(p, lu p ) 
and C(p,u p ) can be extracted from Eq. (0 by suitable 
projections in Dirac-space. 

In Eq. © the Casimir factor C F = (N 2 - l)/(2N c ) 
stems from the color trace. For N c = 2 these equations 
have been solved numerically in Refs. [T3iri4] in a trunca- 
tion scheme which used lattice results from Ref. for 
the temperature dependent gluon propagator as input. As 



discussed in the previous section we now have updated 
and refined results for the SU(2) and SU(3) gluon propa- 
gator at our disposal. These will be used in the following 
to provide for improved results for the quark propagator 
and the associated conventional and dual condensates. 

In the Dyson-Schwinger equation ([6]) for the quark 
propagator we need to evaluate the gluon propagator at 
values of momentum which are not identical with the ones 
extracted on the lattice. In particular we need to evaluate 
the gluon at momenta larger than available lattice mo- 
menta. As it turns out, it is possible to use our knowledge 
on the analytical form of the gluon dressing functions Zt.l 
at zero temperature |43) to devise a fit function which is 
capable to nicely interpolate and extrapolate the lattice 
data also at finite temperature. These fits for the trans- 
verse and longitudinal dressing functions Zt,l of the gluon 
propagator are given by 



ZT,L(q,u q ,T) 



q 2 A 2 



l_ / /?oa(/i)ln[ g 2 /yl 2 + l]^ 
A 2 \ 4tt 



(7) 



Here we find a temperature independent scale parameter 
A =1.4 GeV and the coefficient c = 11.5 GcV 2 . Further- 
more Po = 117V c /3 and 7 = —13/22 in the quenched the- 
ory and we renormalize at = 0.3. 

The fit function Eq. generalises the one used in 
Refs. P~5TfT4^ . where the temperature dependent exponent 
kr,L(T) has been kept fixed, i.e. bT,h{T) = 2. Since we 
now have more accurate lattice gluon data at our disposal 
we found it useful to relax this condition and thus provide 
an even better representation of the lattice data. One may 
speculate whether an irrational exponent bx^hiT) with the 
corresponding temperature dependent cut in the p 2 -plane 
signals quantitative changes also in the analytic structure 
of the gluon, as suggested in [2~4"1|2"5] . A detailed inves- 
tigation of this aspect is relegated to further systematic 
studies [iO] . 

In general our fits are optimised in particular in the 
mid-momentum regime, which is most important later on 
when we use them as input into the quark-DSE. The de- 
tails for the fit parameters as well as plots of the fits for 
selected temperatures are relegated to appendix [D] Here 
we only show the resulting electric and magnetic screen- 
ing masses of the gluon as extracted from the infrared 
behaviour of our fit functions. These are compared to the 
corresponding lattice results in fig. |H The fit quality is 
very satisfactory, except for the electric screening mass 
at the largest temperatures. Here the lattice results suf- 
fer mostly from systematic artifacts due to the restricted 
number of points in the time direction [TD]. This may be 
reflected in the mismatch of the screening masses with the 
ones extracted from our fit. Indeed, one may even argue 
that the electric screening masses from the fits are more 
accurate in the large temperature regime since they nicely 
reproduce the expected proportionality of the screening 
mass with temperature, <~ T, known from hard ther- 
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Fig. 4. Magnetic (left) and electric (right) temperature dependent scale parameters in our fits for the lattice gluon propagator 
compared with the lattice results for the magnetic and electric screening masses of the gluons. 



mal loop results. We also note that at the two highest tem- 
peratures available the fit function describes the low and 
mid momentum behavior of the electric gluon propagator 
very precise. On the other hand the fit function Eq. ([7]) 
is not capable to describe the qualitative mid-momentum 
dependence of the magnetic propagator in this tempera- 
ture range, as can be seen in fig. [9] in the appendix. We 
checked that this momentum behavior can be described 
precisely using in addition a momentum dependent screen- 
ing term in the fit function. Anyway, in the important 
region around the critical temperature both fits work per- 
fectly well and represent therefore a trustablc input for 
the Dyson-Schwingcr equation of the quark propagator. 

Note that as a significant difference to the fits used in 
Refs. [T31H4] it turns out that the transition of the elec- 
trical screening mass from its decreasing behaviour be- 
low the critical temperature T c to the increase above T c 
is much sharper than the one extracted in Refs. |13lH~4] . 
This sharp change around T c was not resolved by the then 
available lattice data of Ref. [T!|| . As a consequence of the 
much improved temperature resolution available now we 
will see that the corresponding deconfinement transition 
extracted from the quark propagator is also much more 
pronounced than the one seen in [T3HT4"1 . This will be dis- 
cussed in more detail in the next section. 

The remaining piece to be specified in the quark-DSE 
is the dressed quark-gluon vertex. Similar to Refs. Q~3, 14 
we employ the following temperature dependent model 



C{k) + C(p) . A(k) + A(p) 
r u {q, k, P ) = Z 3 ( (54,74 - + <E 3 v7j o 



di 



d 2 + q 2 A 2 + q 2 
' /3oa(fi)\n[q 2 /A 2 + l} 

47T 



(8) 



where q = (q,oj q ) denotes the gluon momentum and p = 
(p,Lo p ), k = (fejWfc) the quark and antiquark momenta, 
respectively. Furthermore 25 = —18/44 is the anoma- 
lous dimension of the vertex. Both together, the gluon 
dressing function and the quark-gluon vertex behave like 
the running coupling at large momenta; this is a neces- 
sary boundary condition for any model interaction in the 
quark DSE. The dependence of the vertex on the quark 
dressing functions A and C is motivated by the Slavnov- 
Taylor identity for the vertex; it represents the first term 
of a generalization of the Ball-Chiu vertex [33] to finite 
temperatures. The remaining fit function is purely phe- 
nomcnological, see e.g. |45| where an elaborate version of 
such an ansatz has been used to describe meson observ- 
ablcs. We use d 2 = 0.5 GeV 2 for both gauge groups, but 
dt = 7.6 GeV 2 for SU(2) and d x = 4.6 GeV 2 for SU(3). 
The change in parameter d\ from SU(2) to SU(3) is again 
motivated by the Slavnov- Taylor identity. At high tem- 
peratures it is expected that it reduces to the QED Ward- 
Takahashi identity multiplied with the non-perturbative 
ghost dressing function. A comparison of SU(2) and SU(3) 
ghost dressing functions in the infrared calculated on the 
lattice shows that for small momenta, G(p) of SU(3) is 
reduced by roughly half compared to SU(2), see fig. [5] 
Even though the quantitative values for the ghost dressing 
functions at low momenta from the lattice might contain 
considerable uncertainties we assume the ratio of SU(2) 
to SU(3) to be reliable. In addition we also checked that 
a moderate variation of these parameters does not shift 
the critical temperatures of both, the chiral and the de- 
confinement transition. 

Finally we wish to repeat a word of caution as con- 
cerns the chiral limit in our approximation scheme [13] . A 
prominent feature of the quenched theory not reproduced 
by our framework is the appearance of quenched chiral log- 
arithms in the chiral condensate. These are well-known to 
be generated by rj hairpin diagrams, which are not repre- 
sented by our vertex ansatz. For the present investigation 
we believe this is more an advantage than a drawback. 
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Quenched chiral logarithms are most notable in the chiral 
limit, where they lead to a singularity in the chiral con- 
densate. Since we do not encounter this singularity we are 
in a position to investigate both, the ordinary condensate 
and the dressed Polyakov loop also in the chiral limit. 

The numerical methods needed to solve the quark-DSE 
have been explained in detail in Ref. [14] . Having specified 
all necessary input we now proceed to present our results 
in the next subsection. 



3.3 Numerical results 

Our numerical results for the ordinary quark condensate, 
Eq. (Ql with <p = 7T, and the dressed Polyakov loop, Eq. ([5]), 
are shown in Fig. [SJ Let us first concentrate on the results 
for the ordinary quark condensate. Both, for SU(2) and 
SU(3) we find chiral transitions taking place on a very 
small temperature interval. This is particularly clear in 
the chiral limit (lower panel) but also for quark masses as 
heavy as the strange quark mass (upper panel) . This tem- 
perature interval clearly coincides with the one identified 
in section 12.21 for the change of behaviour in the electric 
and magnetic screening masses. Technically what happens 
in the quark-DSE is that below T c the electric, longitu- 
dinal part of the gluon propagator increases dramatically 
and therefore provides for extra interaction strength in the 
quark DSE. As a consequence, we find increasing values 
for the quark condensate. At or below T c the electric part 
of the gluon propagator reaches its maximum and then 
drops sharply around and above T c . In the quark-DSE this 
sudden loss of interaction strength is responsible for the 
dramatic decrease in the chiral condensate. It is satisfying 
to note that a similar behaviour has been observed from 
calculations of the quenched condensate via Cashcr-Banks 
relations on the lattice [46] . This agreement gives us confi- 
dence that the temperature dependence of our truncation 
for the quark-gluon vertex, Eq. © , is at least qualitatively 
reliable and leads to meaningful results. 

In the chiral limit we clearly obtain a chiral phase tran- 
sition from the conventional quark condensate. Unfortu- 
nately, the temperature resolution of the lattice input is 
still not fine enough to unambiguously identify the or- 
der of the phase transition. One may speculate whether 
the SU(2) transition is second or first order, whereas the 
SU(3) one seems to be first order. The behaviour at finite 
quark masses may be compatible with a rapid cross-over 
for SU(2) and in the case of SU(3) even with a jump in the 
condensate signaling a remnant of a first order transition. 
Further investigations are necessary to clarify, whether the 
differences seen in Fig. [5] between SU(2) and SU(3) are in- 
deed significant. 

As concerns the dressed Polyakov- loop we clearly find a 
transition between the center-symmetric low temperature 
phase and the center-broken phase at transition temper- 
atures very close to the ones encountered for the conven- 
tional quark condensate. Below T c the dressed Polyakov- 
loop is almost constant and very small. For large quark 
masses close to the transition temperature we even find 
small negative values of the Polyakov-loop. We interpret 



these as artifacts introduced due to mass dependencies 
in the quark-gluon vertex that are not represented by 
our vertex ansatz. At temperatures T c < T < 1.1T C the 
Polyakov-loop rises sharply and then less steeply for larger 
temperatures. Within the temperature range investigated 
we do not yet see a saturation of the dressed Polyakov- 
loop at large temperatures, although the results in the 
chiral limit may bear some signals of such a behaviour. In 
general, the deconfinement transition extracted from the 
dressed Polyakov-loop is as pronounced as the correspond- 
ing signal in the electric and magnetic screening masses of 
the gluon propagator, discussed in section HOl Similarly to 
the chiral transitions, also here the temperature resolution 
and the systematic quality of the input gluon propagator 
is not yet good enough to cleanly identify an order of the 
deconfinement transition. 

For SU(2), the results shown here replace previous ones 
reported in Refs. pL"3"l[T4"] . As already discussed above, the 
lattice data used as input in these works have been avail- 
able only for four different temperatures and had to be in- 
terpolated in between. This procedure generated a smooth 
behaviour of the gluon propagator around the critical tem- 
perature resulting also in a broad transition range as con- 
cerns the chiral condensate. Here, with our much better 
temperature resolution as concerns the lattice gluon, we 
are able to identify this behaviour as an artifact of the 
interpolation procedure. 

4 Summary 

In this work we combined two non-perturbativc methods, 
lattice calculations and functional methods using Dyson- 
Schwinger equations, in a common framework to play on 
their individual strengths and to reduce the inherent prob- 
lems of each approach. In the lattice ab initio framework 
we determined the temperature behaviour of the electric 
and magnetic parts of the Landau gauge gluon propagator 
with fine enough temperature resolution to relate their be- 
haviour to the critical temperatures of the deconfinement 
transition of Yang-Mills theory. For SU(2) and SU(3) we 
found a clear signal of the phase transition in the extracted 
electric screening mass and a less pronounced indication in 
the magnetic screening masses. We also verified the earlier 
expectations [TH] that the 'bump' in the electric dressing 
function undergoes dramatic changes around T c . These 
changes have been identified as the source of equally dra- 
matic changes in the ordinary quark condensate as deter- 
mined from the Dyson-Schwinger equation for the quark 
propagator: the condensate keeps rising below T c only to 
decrease sharply at T c . This behaviour is in agreement 
with previous observations from lattice calculations of the 
Cashcr-Banks relation [IB]. It also locates the mechanism 
for chiral symmetry restoration in Landau gauge rather 
unambigously in the (ultra-) soft electric sector of the the- 
ory, emphasizing its genuine non-perturbative nature. 

As for the deconfinement transition we observe a clean 
transition from the dressed Polyakov-loop as determined 
from the Dyson-Schwinger equation for the quark prop- 
agator. The corresponding transition temperatures, T c rj 
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SU(2) massive case^ 



SU(3) massive case] 





T/T„ 



T/T„ 




T/T, 



SU(3) chiral limit 

0.15 




T/T, 



Fig. 5. The quark condensate (tptp)^ and the dressed Polyakov loop Si as a function of temperature for SU(2) (left panel) and 
SU(3) (right panel) Yang-Mills theory. Shown are results in arbitrary units (a.u.) for a massive (strange-)quark (upper panel) 
with m ~ 80 MeV and in the chiral limit (lower panel). 



303 MeV for SU(2) and T c « 277 MeV for SU(3) agree 
with the ones extracted from our lattice calculations. For 
SU(2) our results for the dressed Polyakov-loop update 
and refine the ones of Refs. [T3"irH] . where only a much 
courser lattice input for the gluon propagator has been 
available. 

Finally we wish to emphasize that the dressed Polyakov- 
loop examined in our work is not the only gauge invari- 
ant order parameter for the deconfinement transition that 
can be calculated with functional methods. In Ref. [47] 
the Polyakov-loop potential has been determined within 
the functional renormalization group and delivered the ex- 
pected first and second order phase transition for SU(3) 
and SU(2), respectively. 

In general we believe that all these results demonstrate 
that combined efforts of the lattice and the functional 
framework provide a sound basis to determine gauge in- 
variant observables in gauge fixed Yang-Mills theory. They 
also provide a good starting point to generalise the ap- 



proach towards the inclusion of unquenching effects and 
the introduction of finite chemical potential. 
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A Simulation parameters and volume effects 

In table [T] we present the simulation parameters for the 
lattice calculation of the temperature dependent gluon 
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Table 1. Data of the lattice gauge theory calculations as explained in the text. 



Group 


T/T c 


T [MeV] 


N t 




P 


a [fm] 


Vs' s [fm] 


conf. 


Sill - ?! 


n 
u 


n 


24 


24 


2.227 


fl 1 Q7 


4.74 


V^P, 
oou 


SU(3) 


o 





18 


18 


5.642 


0.197 


3.54 


142 


SU(2) 


o 


o 


24 


24 


2.301 


0.162 


3.89 


189 


SU(3) 


o 


o 


18 


18 


5.738 


0.162 


2.92 


82 


SU(2) 


0.361 


109 


10 


32 


2.261 


0.181 


5.78 


169 


SU(3) 


0.361 


100 


10 


24 


5.642 


0.197 


4.73 


94 


SU(2) 


0.440 


133 


10 


32 


2.332 


0.148 


4.75 


194 


SU(3) 


0.440 


122 


10 


24 


5.738 


0.162 


3.89 


73 


SU(2) 


0.451 


136 


8 


36 


2.261 


0.181 


6.51 


115 


SU(3) 


0.451 


125 


8 


26 


5.642 


0.197 


5.12 


73 


SU(2) 


0.549 


166 


8 


36 


2.3315 


0.149 


5.35 


152 


SU(3) 


0.549 


152 


8 


26 


5.738 


0.162 


5.12 


79 


SU(2) 


0.603 


182 


6 


40 


2.262 


0.180 


7.21 


140 


SU(3) 


0.603 


167 


6 


30 


5.642 


0.197 


5.91 


62 


SU(2) 


0.733 


222 


g 


40 


2.332 


0.148 


5.93 


326 


SU(3) 


0.733 


203 


6 


30 


5.738 


0.162 


4.86 


64 


SU(2) 


0.903 


273 


4 


46 


2.2615 


0.181 


8.30 


246 


SU(3) 


0.903 


250 


4 


34 


5.642 


0.197 


6.70 


59 


SU(2) 


0.968 


293 


4 


46 


2.2872 


0.168 


7.75 


256 


SU(3) 


0.968 


268 


4 


34 


5.675 


0.184 


6.16 


63 


SU(2) 


0.986 


298 


4 


46 


2 2938 


0.165 


7.61 


235 


SU(3) 


0.986 


273 


4 


34 


5.685 


0.180 


6.12 


66 


SU(2) 


1.00 


303 


4 


46 


2.299 


0.163 


7.50 


292 


SU(3) 


1.00 


277 


4 


34 


5.69236 


0.178 


6.05 


65 


SU(2) 


1.02 


308 


4 


46 


2.306 


0.160 


7.35 


280 


SU(3) 


1.02 


282 


4 


34 


5.7 


0.175 


5.95 


67 


SU(2) 


1.04 


314 


4 


46 


2.313 


0.157 


7.21 


319 


SU(3) 


1.04 


288 


4 


34 


5.71 


0.171 


5.81 


113 


SU(2) 


1.10 


333 


4 


46 


2.333 


0.148 


6.80 


276 


SU(3) 


1.10 


305 


4 


34 


5.738 


0.162 


5.51 


71 


SU(2) 


1.81 


548 


2 


52 


2.263 


0.180 


9.35 


232 


SU(3) 


1.81 


500 


2 


38 


5.642 


0.197 


7.49 


59 


SU(2) 


2.20 


665 


2 


52 


2.332 


0.148 


7.71 


310 


SU(3) 


2.20 


609 


2 


38 


5.738 


0.162 


6.16 


70 



propagator. The scale for SU(2) and SU(3) has been fixed 
according to [37] and [35]: respectively, setting the string 
tension to a — (440McV) 2 , to which the considerations of 
[35] apply. T c is then 303 MeV for SU(2) and 277 MeV 
for SU(3), at (3 = 2.299 and j3 = 5.69236, respectively, for 
the employed value of Nt — 4 at T c . T is the correspond- 
ing temperature, Nt the temporal extent of the lattice, 
N s is the spatial extent, (3 the bare coupling, a the lattice 
spacing and V s the corresponding spatial physical volume. 
The number of configurations is denoted by conf., and the 
number of thermalization sweeps is given by 200 + 10iV s , 
and of decorrelation sweeps by 20+ N s . Use has been made 
of the fact that the larger the number of generators, the 
less statistics is needed for the gluon propagator, due to 
the generator-averaging. The statistics was aimed at be- 
low the ten-percent la statistical error level for the electric 
screening mass. 

An immanent drawback for the comparison of SU(2) to 
SU(3) is the lower phase transition temperature of SU(3) 
in physical units. As a consequence, the already more ex- 
pensive calculation for SU(3) would become even more ex- 



pensive if the same physical volumes in spatial direction 
at a fixed time-extension would be required for SU(3). In 
this first investigation for SU(3), this has not been done, 
at the expense of larger finite- volume artifacts at low mo- 
menta for SU(3). These are exemplified in figure [6j where 
the results for the zero-temperature propagator for the 
smallest and largest discretization for SU(3) and SU(2) are 
compared, as an estimate for the systematic uncertainties 
involved. At finite temperature, the spatial volumes are 
larger, but in particular for the chromoelectric gluon this 
is counteracted by discretization artifacts [T9] , 



B Soft vs. hard modes 

At the temperatures studied here, functional studies sug- 
gest that the higher Matsubara frequencies can be ap- 
proximated rather well by Z>r(Wp,p 2 ) = Dt{0,^p +p 2 ), 
and likewise for the other propagators }19il24j . That this 
is indeed a rather good approximation for the volumes 
employed here is shown in figure [JJ 
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SU(2) Gluon propagator at T=0 MeV | 
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SU(3) Gluon propagator at T=0 MeV | 
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Fig. 6. The gluon propagator at zero temperature. SU(2) is shown in the left panel and SU(3) in the right panel. 



Transverse gluon propagator for SU(2r| 
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Transverse gluon propagator tor SU(3^| 
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Fig. 7. Comparison of the soft mode to the first two higher Matsubara frequencies at T = T c . In the upper panels the transverse 
gluon propagator is shown, and in the lower panel the longitudinal ones. In the left panels the results for SU(2) are compared 
to the results for SU(3) in the right panels. All results are shown as a function of the four momentum p 2 = ui 2 + (p) 2 . 



C Ghost dressing function [5] If the ghost is indeed dominated by the gauge-fixing 

procedure this shows that the procedure is not affected 
The soft mode of the ghost is also associated with the s ig nincan tly by a thermodynamic environment. The only 
chromomagnetic sector of the theory [Il[24j[25]. Conse- observable effect, an increase in infrared strength with in- 
quently, in accordance with the expectation [jjj] , it shows 
essentially no dependence on the temperature, see figure 
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SU(2) ghost dressing function below T 
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SU(3) ghost dressing function below T 
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Fig. 8. The ghost dressing function G(p) = p 2 Dc(p) of the ghost propagator Da at temperatures below (top panel) and above 
(middle panel) the phase transition. The bottom panel shows both results together. On the left results for SU(2) and on the 
right results for SU(3) are shown. Lines are drawn to guide the eye. 



creasing temperature, is likely only an effect due to the D Fit functions for the 

larger spatial volumes used at higher temperatures. temperature-dependent gluon propagator 



As already detailed in the main body of this work around 
Eq. Q we use a fit function for the temperature dependent 
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Table 2. Temperature dependent fit parameters for the SU(2) and SU(3) gluon propagator. 
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T/T c 





0.361 


0.44 


0.451 


0.549 


0.603 


0.733 


0.903 


0.968 


0.986 


1 


1.02 


1.04 


1.1 


1.81 


2.2 


a L (T) 
b L (T) 


1.22 
1.94 


0.62 
1.48 


0.48 
1.41 


0.63 
1.52 


0.30 
1.31 


0.29 
1.29 


0.21 
1.29 


0.15 
1.28 


0.15 
1.30 


0.15 
1.29 


0.18 
1.28 


0.34 
1.31 


0.56 
1.41 


0.93 
1.56 


2.69 
1.15 


5.00 
1.16 


a T (L) 
br(L) 


1.22 
1.94 


1.31 
1.98 


1.31 
2.02 


1.39 
2.03 


1.31 
1.98 


1.32 
1.92 


1.32 
1.89 


1.33 
1.76 


1.32 
1.75 


1.29 
1.72 


1.24 
1.71 


1.26 
1.76 


1.30 
1.79 


1.37 
1.83 


1.33 
1.48 


1.35 
1.32 
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1.1 
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a L {T) 
bUT) 


0.60 
1.36 


0.42 
1.23 


0.23 
1.14 


0.33 
1.20 


0.19 
1.13 


0.17 
1.08 


0.11 
1.10 
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1.13 


0.082 
1.14 


0.079 
1.14 


0.16 
1.05 
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1.47 
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0.71 
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0.78 
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1.04 
1.60 


1.05 
1.60 


1.67 
1.91 


1.57 
1.81 


1.06 
1.45 


0.54 
1.13 


0.55 
1.14 


0.57 
1.17 


0.63 
1.19 


1.47 
1.49 


1.42 
1.30 
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Fig. 9. Magnetic (left) and electric (right) temperature dependent gluon dressing function from our lattice calculations 
compared with our fits. 



longitudinal and transverse parts of the gluon propagator 
to represent the corresponding lattice data. 

The resulting fit parameters arx{T) and br,L(T) for 
the gauge groups SU(2) and SU(3) are given in tables [2J 
The rather good quality of the fits, given the statistical 
accuracy and systematic errors of the lattice calculations, 
can be seen in figures 0] and |H1 
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